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Abstract 

The model of Dvali, Gabadadze, and Porrati (DGP) gives a simple geomet- 
rical setup in which gravity becomes 5-dimensional at distances larger than 
a length scale Adgp- We show that this theory has strong interactions at 
a length scale A 3 ~ (A^Qp/Mp) 1 / 3 . If Adgp is of order the Hubble length, 
then the theory loses predictivity at distances shorter than A3 ~ 1000 km. 
The strong interaction can be viewed as arising from a longitudinal 'eaten 
Goldstone' mode that gets a small kinetic term only from mixing with 
transverse graviton polarizations, analogous to the case of massive gravity. 
We also present a negative-energy classical solution, which can be avoided 
by cutting off the theory at the same scale scale A3. Finally, we examine 
the dynamics of the longitudinal Goldstone mode when the background 
geometry is curved. 



*On leave from INFN, Pisa, Italy. 



1 Introduction 



The DGP model is the first ghost-free example of a mechanism in which gravity 
can be localized on a 4D brane in a space of infinite transverse volume. It describes 
a theory where 4D general covariance is unbroken, yet the graviton is a metastable 
state. Its main property is that, on the 4D brane, gravity looks 4D at short distance, 
while it weakens at large distance. This property suggest an interesting alternative 
to the standard description of our present-day accelerating universe. In DGP, the 
cosmic acceleration could be due to gravity becoming weaker at large (horizon-size) 
distance, rather than to a positive cosmological constant. Explicit realizations of this 
scenario have been proposed, for instance in Ref. [2]. 

The model can be described by the action 



S DGP = 2Ml I d 5 x V^G R(G) 

JM 

+ / d 4 x \-4M*Kh) + 2M%Rh) 

JdM L " 



XD 



where Ai is a 5D manifold with boundary dAi, G is the 5D metric, 7 is the 4D 
induced metric on the boundary, and K is the extrinsic curvature. 1 In this model 
gravity on the brane looks 4D at distances shorter than 

a - Ml n 2^ 

Adgp - W ,. (1.2) 

For M 4 3> M 5 this can be a macroscopic length, for example the size of the present 
horizon. 

The DGP model is closely related to massive gravity. In fact, the brane-to-brane 
graviton propagator can be written 

D^ P l(P) = D ™T(P, IpIMdop), (1.3) 

where D™ s £™ e (p,m 2 ) is the propagator for 4D massive gravity. The DGP therefore 
shares with massive gravity the 'vDVZ discontinuity' at distances smaller than 
Adgp; the model reduces not to general relativity, but to a scalar-tensor theory of 
gravity, where the scalar couples with gravitational strength (it does not decouple in 
the limit m — > 0). Refs. [5] showed that in the one graviton exchange approximation, 
massive gravity predicts unacceptable deviations in the predicted bending of light by 
the sun. 



1 The boundary can also be treated as an orbifold fixed point. Wc will discuss the relation between 
these approaches below. 
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However, as shown by Vainshtein |3] for massive gravity, the situation is actually 
more subtle. Near a heavy source, the one graviton exchange approximation breaks 
down at very large distances, and he argued that at smaller distances the resumma- 
tion of nonlinear effects restores agreement with general relativity. In DGP, the one 



Rs is the Schwarzschild radius of the source. At smaller distances, it was explic- 
itly shown that that the full nonlinear solution approaches that of general relativity 
El • Because the scale R* is very large for astrophysical sources, it appeas that 
the DGP model may describe our universe [2J. 

The fact that the one-particle exchange breaks down at a distance so much larger 
than i?s suggests that DGP has hidden strong interaction scales. For a large classical 
source with Rs 3> I/M4 the non-linearities at the scale R* can certainly be associated 
to classical physics. On the other hand, for a source with Rs ~ I/M4, corresponding 
heuristically to one quantum of gravitational charge, we expect any non-linearity to 
be due to quantum physics. Based on this qualitative argument we expect strong 
quantum effects to become important at a length scale 



In this paper we show that this is precisely what happens. For Adgp of order the 
Horizon size, A3 ~ 1000 km. At distances shorter than A3, new interactions become 
important, and there seems to be no reason that the theory should agree with general 
relativity. 

An analogous strong interaction is also present in massive gravity. The strong 
interactions can be made manifest using the Stiickelberg trick of nonlinearly realizing 
the gauge invariance broken by the mass term [S]. In massive gravity, the Stiickelberg 
(or Goldstone) fields have strong self-interactions that necessitate a cutoff that goes 
to zero as the graviton mass goes to zero. An analogous phenomenon is familiar for 
massive non-Abelian gauge fields, where the Stiickelberg sector is a non-linear sigma 
model that is strongly interacting at a scale m/g, where m is the gauge boson mass 
and g is the gauge coupling. In the case of massive gravity, Ref. [S] showed that the 
theory becomes strongly interacting at a scale A5 ~ (m Mp) 1 / 5 , or A 3 ~ (m 2 Mp) 1 ' 3 
if the leading strong terms are tuned to be small. Substituting the 'running mass' 
|p|/Adgp into the A 5 cutoff for massive gravity and solving for p also suggests that 
the DGP model has strong interactions at the scale Eq. (jl.4j) . 

In this paper, we study the DGP model in detail to rigorously establish the exis- 
tence of the strong interactions and understand their origin. Following the logic of the 
Stiickelberg trick, we introduce extra pure gauge degrees of freedom to parameterize 



graviton exchange approximation breaks down at distances R* ~ {Rs^z 




/ 3 , where 




(1.4) 
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the strong interactions. We do this by formulating the theory on a space with bound- 
ary, with no a priori boundary conditions on the fields. This reduces to conventional 
orbifold boundary conditions in a particular gauge, but a different gauge choice is 
useful to make the strong interactions manifest. Since DGP is a generally covariant 
theory, it is not surprising that we find that the Stiickelberg mode has a geometrical 
interpretation: it is a 'brane-bending' mode that keeps the induced boundary metric 
fixed. 

We also find evidence for strong interactions at the scale Eq. (jl.4[) at the classical 
level. We show that the DGP model has classical solutions with negative 5D energy, 
with a boundary stress tensor obeying the dominant energy condition. These solutions 
are at the edge of the regime of validity of the effective theory with short-distance 
cutoff given by Eq. (jl.4|) . giving another indication of new physics at that scale. 

We then consider the behavior of the DGP model in the presence of curvature. 
We show that for the case of a positive curvature boundary (de Sitter sign), the 
self-interactions become stronger, and the Goldstone mode becomes a ghost for suf- 
ficiently large curvature. Closely related results have been found for massive gravity 
in Refs. [H]. We also consider the Randall- Sundrum model ^0] with a DGP kinetic 
term. For a DGP kinetic term on the Planck brane, we find no strong interactions, in 
agreement with expectations based on holography. For a DGP kinetic term on the IR 
brane, the radion becomes a ghost if Adgp becomes larger than the 5D AdS length. 

This paper is organized as follows. Section 2 explains the boundary effective 
action we will use as a tool for the case of a toy scalar model. Section 3 uses this 
formalism to compute the boundary action for the DGP model. We find an effective 
action for the Stiickelberg mode and explicitly compute the cubic interactions. We 
discuss the power counting and derive a non-renormalization theorem for the cubic 
interaction. Section 4 describes a negative energy solution. Section 5 extends the 
analysis of section 3 to the case of backgrounds with nonzero curvature. 

2 Boundary Effective Action 

In this section, we describe the formalism we use to obtain an effective action for 
the boundary field in a theory such as DGP. Consider a field theory on a space with 
boundary, and suppose that we are interested in the correlation functions of sources 
on the boundary. We do not impose any a priori boundary conditions on the fields. 
In the path integral, we integrate over arbitrary boundary values of the bulk fields 
weighted by their action. 

It is useful to separate the fields into bulk fields $ and boundary fields 0. Locality 
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of the action means that the path integral can be written as 

Z = Jd[§] d[<f>] e <W*]+W4>]). (2.1) 

Since the only sources are on the boundary we can integrate out the bulk fields 
to obtain a (nonlocal) effective action for the boundary fields. We must therefore 
integrate over all $ with boundary condition 

$|=0, (2.2) 

where '|' indicates evaluation at the boundary. We perform the $ integral semi- 
classically, by expanding about a solution $ to the bulk equations of motion, with 
5$ | = because of the boundary condition. 2 The semi-classical expression for the 
path integral is then 

Z = Jd[(t)}e^ ShA ^ ]+m ), (2.3) 
where the effective action from integrating out the bulk is 



bulk 



£$2 



$' + ••• , (2.4) 



where the path integral over $' = $ — <£> is performed over fields with boundary 
condition $'| = 0. 



£.i Scalar Field Theory 



Let us do a simple example: free massless scalar field theory in a 5D space with 4D 
boundary at y — 0. The action is 3 



(fx 



-\d M <$>d M § 



y=0 



where 



$| = 0. 

The classical solution for $ with these boundary condition is 



-yA 



(f)(x), 



(2.5) 



(2.6) 



(2.7) 



2 Because S$>\ = there is no boundary term in the variation of the bulk action. 
3 If we had written the bulk kinetic term as ^$IZl5<J> there would be a boundary term in the 
variation proportional to <J><9j,(<5<£>), affecting the behavior of solutions near the boundary. 
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where A = y/— O4. We therefore obtain 



1*] = /. 



y=0 



±0A0 

From this we can read off the propagator for the field: 

f 



(00) 



Now we add bulk interactions: 



kUa - A' 



-±<9 M $<9 M - 



1 \^f\3 



Sbulk — Id X 

The classical bulk field satisfies 

□ 5 $_A$ 2 = 0, $| = 0. 
We find the solution order by order in A: 

<5 = <l + $1 H , 

where <J n = 0(X n ). $0 was computed above. Because $ | = 0, we have 

$ n | = for n > 1. 
We now compute the first-order correction to the action: 

r[0] = r + r 1 + ---, 

where T was computed above, and 

r\ = J (fx [-d M ^d M % - k x ® 3 



d 5 x 



3^0 

d 5 x $, 3 



y=0 



TJ A j it ■ I J 1 1 

_i A f d% d 4 p 3 
~ 3 J (2tt) 4 '" (2tt)' 



(2.8) 
(2.9) 

(2.10) 

(2.H) 
(2.12) 
(2-13) 
(2.14) 
(2.15) 



w'av+p.+fto (2.16) 



where 



0(p) = Jd 4 xe tp - X <j)(x). 



(2.17) 
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3 Boundary Effective Action for Gravity 

We now consider 5D gravity with a 4D boundary, the case of interest. We use Latin 
capitals M,N, . . . = 0, . . . , 3; 5 for 5D spacetime indices, and u, . . . for 4D ones. We 
denote the bulk metric by Gmn- 

In the standard treatment, we define boundary conditions by imposing an orbifold 
projection under reflections about the boundary. Here instead we will not impose any 
boundary conditions on Gmn- This means that there is extra gauge freedom in this 
formulation. In the bulk, we have infinitesimal gauge transformations generated by 

SGmn = ^ P 9pGmn + 9m^ f Gpn + 8n^ p Gmp- (3.1) 

In the orbifold formulation, we have G^\ = and hence = 0, H 5 | = 0, where 
the dot denotes the derivative with respect to x 5 and the vertical stroke denotes 
evaluation at the boundary. In the present formulation, G^^\ ^ and 5 M | ^ 0. We 
still have S 5 | = because we use coordinates where the boundary position is fixed. 
We can choose a gauge where G^l = to recover the orbifold boundary conditions, 
so this formulation is completely equivalent to the usual one. However, the extra 
gauge degrees of freedom in the present approach are very useful in uncovering the 
strong interactions, as for massive gravity jH]. 

It is convenient to make a 4 + 1 split and write the action in terms of ADM-like 
variables [UJ: the lapse N = (G 55 ) -1 / 2 , the shift = G 5fJi , and the 4D metric 
Ifiu = G p on surfaces of constant y = x 5 : 

S hulk = 2Ml Jd 4 x 

where 

K fU ,= 

is the extrinsic curvature. Here, 4D indices are raised and lowered with 7^, is 
the covariant derivative with respect to the 4D metric 7^, and the dot denotes a 
derivative with respect to y. Note that only first derivatives appear in the action 
and there is no boundary (Gibbons-Hawking) term in this formulation [TT) IT!?] (see 
also [TH]). 

In order to integrate out the bulk fields we must choose a gauge for them. We 
want to choose a gauge such that the propagator has manifestly good high-energy 
behavior, so we choose de Donder gauge. We write 

Gmn = Vmn + H M n, l^v — Vnv + (3-4) 



dy v^7iV R{ 7 ) - K^K^ + K 1 



(3.2) 



IN 



( V - D„N V - D V N^) 



(3.3) 
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and add the gauge fixing term 

£buik, gf = -MlF M F M , (3.5) 

where 

F M = d N H MN — \d M H. (3.6) 

Classically, this imposes the gauge Fm = 0. In terms of the 4+1 split, 

£buik, gf = -Mf[(c^ - \duC ~ \d v H^ + N v f + (d»N, + \H„ - \Q} 2 (3.7) 

This leaves residual gauge freedom under infinitesimal transformations satisfying 

□ 5 S M = 0. (3.8) 

Because we require S M to be well-behaved at infinity, and S 5 | =0, we see that S 5 is 
completely fixed but there is a residual gauge freedom, parameterized by ^ = S M |. 
Explicitly, the residual gauge freedom is 

3„ = e~y%. (3.9) 

This residual gauge freedom acts on the boundary fields at linear order as 

5h lu/ = d^ v + d v ^, 5N II = -A^, 5h 55 = 0. (3.10) 

We now integrate out the bulk fields to obtain the quadratic boundary action. We 
solve the bulk equations of motion with boundary conditions 

Hmn = hMN- (3-11) 

The equations of motion in de Donder gauge are 

n n (H M N - IvmnH) = 0, (3.12) 

with solution 

H M n = e~ yA h M N- (3.13) 

The induced boundary action is 

T = Ml |d 4 x[-|FAV + i/i 4 A/i 4 + ±h 4 Ah 55 - ±h 55 Ah 55 

- N^AN^ - N»(d^h 4 + d„h 55 - 29" V)l • (3- 14 ) 
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This is invariant under the gauge transformations Eq. 1)3. 10|l . In fact, in terms of the 
invariant combination 



1 



we have 



T = Ml ld*x 



-fh^Ah^ + \h A Ah A + fh A Ah 55 - lh 55 Ah 55 



(3.15) 



(3.16) 



The induced boundary action must be added to the DGP kinetic term on the 
boundary: 



£bdy,DGP — 



(3.17) 



We fix the remaining gauge freedom parameterized by £^ by adding a gauge fixing 
term 



C 



bdy,gf 



-Ml(d»h^-\d v K + mN v ) 



(3.18) 



where m = M^/Mf is the DGP scale. This gauge fixing makes the large DGP kinetic 
term invertible, and also eliminates the mixing between h^ u and iV M . The complete 
quadratic boundary Lagrangian is then 



C 



bdy 



Ml 



|F(n4-mA)V - \h A (U A -mA)h A 



-mN^(A + m)N tl + \mh A Ah bb - mN^d^h^ - m\h 55 Ah 55 
To see the strongly interacting mode, we consider the scalar modes 

h 

and h 55 . In the regime p ^> m the leading terms are 



bdy 



Ma 



>U(f) + 2m<f)Ah 55 - m(a + ~/i 55 )A(cr + \h 



55 / 



(3.19) 



(3.20) 



(3.21) 



From this we see that there is one scalar mode that gets a kinetic term only through 
mixing with h^ u . This mode can be parameterized by 



= <9 M 7T, h 55 = -2A7T. 

We can diagonalize the full kinetic term by defining 



iv; = N, - «Vr, 



ti. 



ft]' 



(3.22) 



(3.23) 
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and we obtain 



£bdy ^ Mi 



^□ 4 /^-^n4/4-mA^AiV; + 3m 2 7rn47r . (3.24) 



The small coefficient of the 7r kinetic term is the origin of the strong interactions in 
this theory. 

We can characterize the strongly-interacting mode in another way, which makes 
the generalization to curved backgrounds more transparent. The mode we found can 
be characterized by the following three properties: (i) it solves the linearized bulk 
equations of motion; (ii) = 0; (Hi) it obeys the de Donder gauge-fixing condition 

F N = d M H MN - \d N H = 0. (3.25) 

To see this, note that and H 55 can be locally gauged away, so any configuration 
satisfying these conditions must be pure gauge in the bulk at linear order. In fact, 
the mode Eq. (|3.22j) extended into the bulk is 

= 0, = dfjE 5 , H 55 = 2S 5 , (3.26) 

where 

E 5 = e~ yA 7r. (3.27) 

Since S 5 | ^ this gauge transformation is not a symmetry of the full action with 
boundary. The boundary shifts under this transformation, so this can be viewed 
as a 'brane bending mode.' This is the only nontrivial configuration with the three 
properties described above. Beyond linear order, S 5 | and 7r have different interactions, 
since S 5 | (unlike 7r) affects the 4D induced metric at order (S 5 |) 2 

In fact, one could have anticipated by purely geometrical and physical arguments 
which mode, if any, could interact strongly. The strong mode should be related with 
the UV properties at the boundary, so we expect it to correspond to a trivial bulk 
geometry away from the brane. This is to say that the mode should be pure gauge 
in the bulk. Moreover, it should also not correspond to sizeable curvature of the 
induced boundary geometry: this is because the large DGP Einstein term disfavors 
intrinsic curvature. Since both the bulk and brane geometry should not be excited, the 
only remaining geometrical object that can be excited by the mode is the extrinsic 
curvature of the boundary, describing its shape as seen by a 5D observer. Up to 
trivial 4D reparametrizations, the only mode satisfying the above three requirements 
is precisely Eq. ()3.26|) . Notice that the second requirement bears similarity to the 
case of massive non-abelian gauge theory. There, the strongly interacting Goldstones 
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are the the pure gauge configurations = Wd^U for which the gauge kinetic term 
vanishes. The gauge kinetic term, whose coefficient 1/g 2 in principle can be very 
large, is the analog of the DGP term. 

We now turn to the question of higher- dimension operators in the effective theory. 
Bulk interactions with higher powers of Hmn will give rise to boundary interactions 
of the form 

A£ bdy ~ Mld(N,r(dnY ~ mM\d [J^] (^] ' , (3.28) 



m 



where 

mM 4 vr, ~ m 1/2 M 4 A^ (3.29) 



7T 



are the fields with unit kinetic term. From this we can read off a strong interaction 
scale 

A M ~ (m^-Hir q - 2 ) 1/(3p/2+2q+1) . (3.30) 

The lowest scale occurs for p = 0, q = 3 (cubic tc interactions), which gives a scale 

A ~ (m 2 M 4 ) 1/3 . (3.31) 

Higher derivative terms in the bulk give rise to terms with additional powers of d/M$. 
Since M5 ^> A, these will give weaker interactions. 

Notice that cubic terms cannot be canceled by changing the gauge condition. 
Suppose that we modify the gauge-fixing condition by adding non-linear terms in 
H MN : F' M = F M + 0(H 2 ). Since the mode Eq. (|3.22|) obeys F M = 0, the new gauge 
condition becomes F' M = 0(n 2 ). This change produces only terms of order 7r 4 or 
higher in the boundary action. 

We now show that the cubic terms are present by computing them explicitly. We 
must evaluate the cubic terms in 7r in the configuration (see Eq. (|3.22jl ) 

= ^n, H 55 = 2«9 5 n, (3.32) 

where 

n = e~ yA 7i. (3.33) 

To find the cubic terms in the bulk action Eq. ()3.2)1 we need A^ to linear order and 
Kiiv to quadratic order: 

N=l + §# 55 (3.34) 
= |(1 - §#55) WV"* + a » N v)- ( 3 - 35 ) 
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The cubic terms involving tt and iV„ are 



A£ bulk = M 5 



^(•V^ - dvNrf + d"H 5S N li ff'N v - d u H^N^N v \ . (3.36) 
Using Eq. ()3.32|) and All = —II, this can be written 

A£ bulk = 2M|(^An)(^nn 4 n - d v m^d v \i) 

= -Mld y [^n^nn 4 n] . (3.37) 

Integrating this solution over y, we obtain 

A£ bdy = Mlmd^Ttd^U^K. (3.38) 

To see that this cubic interaction represents a physical effect, we can compute the 
correlation function of three stress-energy tensors on the brane. For kinematics where 
all momenta are space-like off-shell with p ^> m, the cubic interaction computed above 
dominates the amplitude, which becomes strong at the scale Eq. (|3.31j) . Similarly, 
by studying the Feyman diagrams, one finds that Eq. ()3.38|) leads to a non-trivial 
4-point scattering amplitude that violates unitarity at the scale A. 

When we include quantum corrections, we expect all operators consistent with 
symmetries to be generated, and we expect an infinite number of terms that get 
strong at the scale A. Indeed, for the subset of logarithmically divergent graphs, we 
must include the associated operators in order to be consistent with unitarity. 4 These 
terms must be localized at the boundary, since the cutoff for bulk interactions far from 
the boundary is M 5 ^> Adgp- These interactions must respect 4D Lorentz invariance 
and must be local when written in terms of the geometrical quantities R^pail) arid 
K^v. In order to zoom in on the strong interactions it is convenient to take a limit 
where A = M|/M 4 is fixed and M 4 , M 5 -> oo, so that also m = M|/M| -> 0. This is 
the analogue of the g —>■ with /„. = my / g fixed limit of massive non-abelian gauge 
theory. In this limit, the geometrical objects reduce to their linearized approximation 

m -X,= ^ + 0(^), (3.39) 
m - lK ^Mf + (^m (3.40) 



4 In other words, logarithmic divergences correspond to the RG evolution, so that the coefficient 
of the corresponding operators cannot be set to zero at all scales. Power divergent effects are not 
calculable and could consistently be set to zero, for example by using dimensional regularization. 
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where n = it/M^m is the canonically normalized field. By these equations we expect 
the terms that get strong at the scale A to have the form 



(3.41) 



Note that the theory becomes strongly coupled whenever the 4D curvature is of order 
m 2 . We will comment further on this point in Section 4 below. 

Note also that the cubic interaction Eq. ()3.38|) is nonlocal when written in terms 
of geometrical quantities (since K ~ d 2 ir and R ~ d 2 7i). This means that loops of 
7r fields should not renormalize this interaction, and that the divergent part of loop 
diagrams involving this interaction should be expressible as a function of d 2 ir. 5 This 
non-renormalization theorem follows simply by integration by parts. Consider any 
1PI diagram with an external line coming from one of the factors of 7r with only one 
derivative. Because the diagram is 1PI, both of the other 7r factors attach to internal 
lines. We then have 

^7r ext <9 /1 7ri nt n 4 7r int = d^n cxt d v 9 M 7r int ^vr int - |^9 p vr int 9 p 7r int , (3.42) 

which is a function of <9 2 7r ext after integration by parts. This gives a nice check of the 
consistency of this framework. 



4 Classical Instabilities 

In this section we study a classical solution to the DGP model in which the stress- 
energy tensor on the brane satisfies the dominant energy condition, yet the brane has 
negative energy from the 5D point of view. When the boundary has the topology of 
R 4 it is difficult to define the 5D energy, which is presumably infinite. We therefore 
look for static solution where the spatial sections of the boundary have topology S* 3 , 
the bulk is 'outside' the S s , and the solution is 0(4) symmetric. The geometry of the 
boundary therefore corresponds to a spatially compact static cosmological solution, 
similar to the Einstein universe. 

By Birkhoff's theorem, the metric outside the boundary is the 5D Schwarzschild 
metric 

ds 2 = -f 2 (r)dt 2 + + r 2 dn 2 3 , f = (4.1) 

5 A full calculation should include loops of bulk fields as well. However, the scaling argument 
shows the leading interactions at the scale A are expressible in terms of the interactions of the tt 
field. 
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where Rs is the Schwarzschild radius. The boundary is at a fixed value of r > in 
these coordinates. Because this solution asymptotes to 5D flat space at infinity, the 
energy (mass) of the solution is well-defined: 

M = 327rM 5 3 i?|. (4.2) 

For _R| < such a solution has negative energy. It is a negative-mass Schwarzschild 
solution with the naked singularity cut out by the boundary. 

The most general form of the stress-energy tensor on the brane compatible with 
the symmetries is 

Too = -ploo, Tij = +Plij, (4-3) 
where p is the energy density and p is the pressure. We impose the equation of state 

p = wp. (4.4) 
We look for solutions satisfying the dominant energy condition, which requires 

P > 0, -1 < w < 1. (4.5) 

The bulk Einstein equations are satisfied by the metric Eq. (j4.1|) . The only addi- 
tional equation that must be satisfied is 

4M 4 2 ^( 7 ) - 4Ml(K^ - lvv K) = T^, (4.6) 

where Q^ u is the Einstein tensor and T^ v is a stress tensor on the boundary. Eq. (|4.6J) 
follows simply from the variation of the full action with free boundary conditions. 
(There is no junction equation in this approach since there is no 'other side' to the 
boundary.) In the metric Eq. (j4.1|) . we can use Eq. ()3.3|) to obtain 

K, v = \fd rl ^ v . (4.7) 

The boundary equation (|4.6J) then gives 

AMl--AMf-l = Pl (4.8) 

±Ml±-4M$l(f+j^=-wp. (4.9) 

Note that when M 5 = the 4D solution reduces to the standard Einstein static 
universe with w = —1/3, p = 12M|/r 2 . When M 4 = 0, Eq. ()4.(ij) is equivalent 
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to the usual Israel junction conditions, and Eqs. (|4.8J) and (|4.9J) have no solutions 
satisfying the dominant energy condition. (In fact, they have no solutions even if the 
stress energy tensor is allowed to be written as a negative tension term plus a term 
satisfying the dominant energy condition.) 

The constraint p > gives 

rf < — . (4.10) 

m 

Combining Eqs. f)4.8|) and ()4.9|) we obtain 

1 4M. 3 

• = -3 + -^. (4.11) 
which shows that w > — 1 is satisfied whenever p > 0. The condition w < 1 gives 

r/ + ^7<-- (4.12) 
4/ m 

Since this is clearly more restrictive than Eq. ([4.10)1 . this is the only condition that 
needs to be checked. For M < (i?| < 0), the left-hand side of Eq. f)4. 12|) approaches 
|i?s| as r — > and increases monotonically with r, so we obtain a solution for \Rg\ < 
1/m. This means that the energy cannot be made arbitrarily negative in this model 
(see Eq. (jOJ). 

Note that the minimal 4D curvature of a negative energy solution is 0(m 2 ). For 
the critical zero energy solution / — > 1, and Eq. (|4.12|) implies r < 1/m, so > m 2 . 
However, Eq. ()3.41|) shows that such curvature is also the critical curvature where 
the derivative expansion of the effective quantum field theory breaks down. So the 
negative energy solutions to lie at the edge of the regime of validity of our theory. 

Another way of arguing the same point is the following. The instability appears 
for 4D energy density p = 0(M%m 2 ). Significantly, this is also the energy density 
of a gravitational source for which the cubic interaction term in Eq. ()3.38j) becomes 
comparable to the kinetic term of the scalar tt, defined in Section 3. To see this, recall 
that 7r couples to the stress-energy tensor with strength m (see Eq. (|3.23|) ). Then, 
to linear order in the source, ~ T^/M 2 m. Substituting this estimate into the 
cubic interaction term Eq. ()3.38|) we have 

A£ bdy ~ T\d v nd v n. (4.13) 

This term becomes of the same order as the 7r kinetic term when T M M ~ M|m 2 . We 
conclude that the negative energy solutions appear only at the edge of validity of 
the effective theory, and the theory with cutoff of order A ~ (rr^M^) 1 ^ 3 is safe from 
instabilities. 
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5 Curved Backgrounds 



A noteworthy aspect of massive gravity is that when propagating on a curved back- 
ground, it behaves very differently than in fiat space. In AdS space, there is no 
vDVZ [3| discontinuity [9. , while in dS a light massive graviton becomes a ghost |14j . 
These unusual features find a simple explanation when massive gravity is rendered 
covariant by adding a Goldstone vector [T51IE] A^. At linear order, A^ appears in the 
combination h^ v — DuA v ) . Here is the covariant derivative of the 4D background. 
The difference with flat space originates from the fact that at nonzero cosmological 
constant A, the kinetic term of the strongly-interacting scalar mode 7r inside A^, 
Afj, = D^tt, receives a contribution proportional to —A. This contribution suppresses 
the cubic interactions of 7r when A < 0, and makes the graviton a ghost at small mass 
when A > 0. 

In DGP, we find an analogous phenomenon. We first give a general argument, then 
consider two important special cases: one with a boundary with de Sitter geometry, 
the other a Randall-Sundrum model with a DGP kinetic term on the boundaries. 



5.1 General Discussion 

We consider the linearized theory about a general 5D background metric 

Gain = Gmn + H MN . (5.1) 

We generalize de Donder gauge to curved backgrounds by adding the gauge fixing 
term 

A£ bulk , gf = — M| G mn Fm Fn , (5.2) 

where 

Fm = V Hmn — |V mH, (5.3) 

where Vm is the covariant derivative associated with the background metric Gmn- 
The bulk equations of motion are then 

V 2 (Hmn - \GmnH) = 0. (5.4) 

Following the discussion in the flat case, it is convenient to parameterize the modes 
of Hmn as follows. For simplicity, we will give the discussion using Gaussian normal 
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coordinates for the background metric: = 0, G55 = 1. Instead of H 5fl \ = h 5fl and 
Hsb\ = h 55 , we use as boundary variables 

= ^, S 5 | = 7T, (5.5) 

where "Em satisfies the bulk equation 

V 2 S A/ = 0. (5.6) 

This implies that the mode H MN = V(mSj\t) satisfies the bulk equations of motion as 
well as the de Donder gauge fixing condition. We then parameterize a general bulk 
fluctuation as 

Hmn = H$ N + V(m5jv) — Hmn, (5.7) 

where 

V 2 #Sv = 0, V 2 #Sv = 0, (5.8) 

with boundary values 

Hj$\=h ia „ ^1=0, #g } | = 0, (5.9) 

and 

=V (M E U) I =D^ U) + 2K^ 5 , H 5fl \=0, H 55 \=0, (5.10) 

where D M is the covariant derivative with respect to the induced background metric 
7 M „. The Hmn term in Eq. (|5.7j) subtracts the contribution of the 'Goldstone' modes 
E M to the fluctuations of the induced boundary metric, which is then simply Hj?}\ = 
h^v, the 'spin 2' mode. This ensures that £ M and 7r as defined above do not appear 
in the large DGP kinetic term. We will see that in terms of these variables, the 
identification of the strong degrees of freedom is more direct. (Indeed, the combination 
N' that diagonalizes the kinetic term in Eq. (J3.23)) in the flat case is precisely the 
'Goldstone' 

The dependence on £ M can be obtained simply by noting that under 

hfj, v 1— > D(nX v ), £ m >-> £ m + A m , 7T 1 > 7r, (5.11) 

the bulk fluctuation Eq. (J5.7)) changes precisely by a residual gauge transformation 
A M , satisfying A M | = A M , A 5 = 0. Eq. (|5.11j) is therefore a symmetry of the quadratic 
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boundary action. Therefore, we can work out the action at £^ = and restore the 
dependence on £^ by the substitution 



h 



(IV 



h 



Dfrtv). (5.12) 
We now compute the boundary action for the modes above. For this we will need 
V 5 iJ$|~-AV, V 5 H M |~-A£ M , V 5 n|~-A7r, (5.13) 
valid for modes with 4D wavelengths smaller than the scale of curvature, where 

A = sf^D 2 . (5-14) 

Eq. (|5.13|) follows from the fact that for small-wavelength fluctuations the curvature 
is irrelevant, so the result must reduce smoothly to the flat case. We will see how 
this arises in an explicit calculation in the next subsection. The contribution to the 
boundary effective action from the bulk Einstein action is then 



AC 



bdy,E 



linearized 



-y=fM 5 3 [irA/i F - \hAh + brDyPv-R - hD\ 



(5.15) 
(5.16) 



where we have omitted terms of order Kh 2 and KhAir that are subleading at small 
4D wavelengths. The contribution from the bulk de Donder gauge fixing term is 



A£bd y ,buikgf = — ^1 v / — 7 \ H 5 mF m — \HF" 



-\hAh 



ir(D^D u h^ - D 2 h) 
2n(K flu D fl D u - KD 2 )n + 



where 



= - 2K^n) - \D^{h - 2Kit) + 0(K 2 tt, Kh), 

F 5 = \A(h - 2Kt() + 0(K\, Kh), 



(5.17) 
(5.18) 

(5.19) 
(5.20) 



and we again omit terms that are subleading at small 4D wavelengths. At this point, 
it is trivial to include the mode ^ by the substitution Eq. ([5. 12)1 . Note that n mixes 
with h^ u only via the linearized Ricci scalar, so there is no mixing between 7r and £ M , 
generalizing the result of flat case. 

The full quadratic boundary action is therefore the sum of Eqs. ()5.16|) and (|5.18|) . 
In the small wavelength approximation, the only relevant change comes from the last 
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two terms in Eq. (j5.18j) . To understand what they imply, we consider for simplicity 
a maximally symmetric background, for which 

K^u = C% v . (5.21) 

Adding the DGP kinetic term and diagonalizing the kinetic term for 7r and we 
obtain the 7r kinetic term 

£ kin = 3M 4 2 m(m - 2C)nD 2 n. (5.22) 

Note that the n mode becomes a ghost for C = |m. We will see below that C > for 
de Sitter space (see Eq. ()5.27|) ). Therefore, as in massive gravity, positive curvature 
increases the strength of the strongest interactions, while negative curvature decreases 
it. 

One of the motivations for considering the DGP model is that it provides a source 
of 'dark energy' in the absence of 4D vacuum energy In this case, we can compute 
the relation between the constant C introduced above and the positive curvature of 
the present-day universe. From Eq. 1)4.6)1 . we have 

4M 2 £^( 7 ) = - 12M 5 3 CV ~ - 4M 2 A now7 ^. (5.23) 

This shows that C is positive, hence 7r is a ghost, for sufficiently small m. In the next 
subsection, we will see that in a DGP model with no cosmological constant in the 
bulk, 7r is a ghost in the regime where 4D vacuum energy does not contribute to the 
4D curvature. More generally, we expect that positive 4D curvature decreases the 
strength of the 7r kinetic term, so that it makes the interactions of 7r even stronger 
than on a Minkowsky 4D background. Conversely, a negative 4D curvature weakens 
the 7r self-interactions, as in massive gravity. When m is much smaller than the 
curvature |C|, it even eliminates the vDVZ discontinuity already at linear order jH]. 

5.2 Explicit Calculation: de Sitter Space 

We now consider the important special case of 4D de Sitter space in a DGP model with 
vanishing bulk cosmological constant. The solution has very simple 5D geometry: the 
bulk is flat, and the boundary is at 

V^xV + y 2 = L 2 (5.24) 

in Cartesian coordinates. It is more convenient to use coordinates where the back- 
ground metric is 

ds 2 = dr 2 - r 2 dr 2 + r 2 cosh 2 r dtt 2 3 . (5.25) 
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The boundary is at r = L in these coordinates. The boundary equation is 

4M 4 2 £ Ml ,(7) - 4M 5 2 [K^ - % U R] = -V %„, (5.26) 

where Vo > is the vacuum energy on the boundary. It is straightforward to work 
out 

3 1 

Gnv(l) = Kp, = (5.27) 

This gives the relation between the 4D vacuum energy and the curvature: 

V 



Y1M\ 



L 2 + mL = l. (5.28) 



For VqI M\ > m 2 , this gives the usual relation between vacuum energy and curvature, 
but for Vo/Mj -C m? we get L ~ 1/m [2]. It is interesting that in this regime the 4D 
curvature is independent of the vacuum energy, but from Eq. ()5.22j) we see that the 
Goldstone is always a ghost in this regime. 

We now consider the explicit calculation of the boundary action. To understand 
the strong interactions, it is sufficient to consider the scalar modes 

H M n = G M n® + (n M VjvS + %VmS) + n M n N Q, (5.29) 
where %/ is the normal vector n 5 = n 5 = 1, n M , = 0. A useful relation is 

Vm^tv = -Imn, (5.30) 
r 

where 

Imn = Gmn ~ n M n N (5.31) 

is the induced metric on surfaces of constant r. 6 The de Donder equations of motion 
for the scalar modes defined above are 

%» (v 2 $ + ^Q^j + * V M V„E = 0, (5.32) 



2 _ 4 9 

v 2 s + -v r E - — e + -n 



ry ry 2 ry 

4 

,y2 J ry£. 



0, (5.33) 



V 2 fi + V 2 $ + 2V r V 2 E- — E -— Q = 0. (5.34) 



6 We can think of r as a bulk scalar. Geometrically, it is the proper distance of any point in the 
bulk to the boundary. 
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In our coordinates 



V 2 $ 



4 1 
d* + -d r + — a 



$, (5.35) 



where CU is the Laplacian on S 3 . Since D 4 is independent of r, we can treat it 
as a parameter when solving the equations. (More formally, we could expand in 
eigenstates of D 4 .) 

The behavior of the solutions is very easy to understand once we notice that they 
are homogeneous in r. The solutions therefore have the form 



where A depends on n 4 . The solution is very simple for fluctuations with D 4 ^> 1/L 2 
on the boundary. (Recall that L is the size of the 4D universe!) For these fluctuations 
□ 4 ^> 1 and the leading terms in the equations that determine A are simply 

A 2 + D 4 = 0, (5.37) 



with solution A = ±y — n 4 . Good behavior at infinity (away from the boundary) 
requires the negative solution. Since n 4 is the only expansion parameter, the correc- 
tions are 



A = -V-D4 [1 + C(1/D 4 )J , (5.38) 

and so on the boundary 

A\ = -LA + 0(1/ LA). (5.39) 

Note that the corrections to A are smaller than the 1/L curvature corrections de- 
scribed in the previous subsection. 

The conclusions depend only on the fact that the equations are second order and 
homogeneous in r, and so hold for general polarization states. This shows that 

d r H MN \ = -AH MN + 0(H/L 3 A). (5.40) 

The remainder of the calculation follows the previous subsection line by line. 

5.3 DGP and Randall- Sundrum 

It is instructive to apply the results of the previous section to the Randall-Sundrum 
(RS) model [TU] with a DGP boundary term added. At the boundaries, we have 

= ±L% U , (5.41) 
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where L — 1/k is the bulk AdS curvature length, and the + (— ) sign corresponds to 
the IR (Planck) brane. 

We first consider adding a DGP kinetic term on the Planck brane. As long as 
k > m, Eq. ()5.22j) shows that the scale of strong interactions is (M 5 /c) 1//2 > k. We 
conclude that the DGP kinetic term does not lead to new strong interactions within 
the original regime of validity of the RS model. This is consistent with the holographic 
interpretation of the model as a 4D conformal field theory (CFT) coupled to gravity. 
The DGP kinetic term simply corresponds to a large coefficient for the gravity kinetic 
term in the UV, large enough to dominate the induced contribution to the Planck 
scale from the CFT. Note that the extrinsic curvature term in the boundary effective 
action is crucial for obtaining this result. 

We now consider adding a DGP kinetic term to the IR (or 'TeV') brane. In 
this case, Eq. 1)5.22)1 tells us that the theory has a ghost for k > |m. In fact, this 
instability comes from the radion itself becoming a ghost. To see this in a simple 
way, we consider the limit where the Planck brane is pushed to the boundary of AdS 
and the IR brane is at fixed position y = 0. In this limit the 4D zero mode graviton 
is decoupled from the physics on the IR brane. The only zero mode is the radion 0, 
which can be parameterized by JE] 



in agreement with Eq. (j5.22j) . 

This is interesting because it prevents a geometrical construction of a model with 
an isolated massive spin 2 particle. This model has only massive spin 2 KK modes, 
and when M| ^> M^/k one finds that the lightest spin 2 mode has an anomalously 
small mass of order (km) 1 / 2 , while the remaining massive KK modes have mass of 
order k. If it were not for the instability, discussed above, for m <C k this model 
would reduce to an effective theory of a single massive graviton (and a radion) below 
the scale k. 



ds 2 = e 2ky \l + 2<j)(x)e- 2ky ] dx^dx 1 " + \l - 2(j)(x)e~ 2ky ] dy 2 . 



(5.42) 



Its kinetic term is 




(5.43) 
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6 Conclusions 



We now summarize our results. First, we showed that the DGP model has strong 
interactions at distances shorter than A3 ~ (-^dgp/^) 1 ^ 3 - The strong interactions 
are due to a scalar 'Goldstone' mode that obtains a kinetic term only by mixing with 
the transverse graviton polarizations, similar to massive gravity. The longitudinal 
Goldstone has a geometrical interpretation as a brane bending mode that keeps the 
induced metric on the brane fixed. Second, we showed that there are classical insta- 
bilities in the DGP model in the form of negative energy solutions. These solutions 
are at the edge of validity of the effective field theory with UV cutoff at the scale A3, 
giving further support to the conclusion that new physics is required at this scale. 
Finally, we considered the strong interactions in the presence of 4D and/or 5D cur- 
vature. We showed that positive (de Sitter) sign curvature makes the model more 
strongly interacting, and makes the strong mode a ghost for sufficiently large curva- 
ture. We also investigated the effect of a DGP kinetic term in the Randall-Sundrum 
model. 

We conclude with some comments on the solution of the vDVZ [3] discontinuity 
problem suggested in Ref. [5J. There, it was shown that around a classical source 
of Schwarzschild radius Rs, a careful resummation of non-linear effects restores the 
phenomenologically correct Schwarzschild solution below a distance R* ~ (Rs/m 2 ) 1 ^ 3 . 
Now, the curvature of the Schwarzschild metric is of order R$/r 3 , so that at the 
distance R* the curvature is of order m 2 . From inspection of Eq. ()3.41|) we find 
that this is the critical curvature at which the quantum effective field theory breaks 
down. Stated otherwise, at the scale R*, the effective quantum expansion parameter 
e = <9 2 7r/A 3 becomes order 1. Therefore, any statement about the behavior of the 
field at distances smaller than R* requires knowledge (or assumptions) about the UV 
completion of the 7r sector. 

One possibility is to assume that when e ^> 1 the infinite series of counterterms 
saturates and the quantum correction to the effective action stays of the order of the 
result at e ~ 1. In this case, the contribution from the counterterms is suppressed 
compared to the tree-level contribution (which has fewer derivatives) by ~ 1/ (M^Rs) 2 . 
In this scenario, there is a range of scales where the effect of Ref. [3] works. In this 
case, one still has to find a way to avoid the instabilities associated with curvature of 
order m 2 , found above. Moreover, this does not rescue the model phenomenologically, 
since at least at the length scale A3 ~ 1000 km, gravity becomes sensitive to the details 
of the UV completion. 
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